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Abstract 



Open, dissipative systems subject to a random force are directly 
quantized. The starting point is the effective action derived using 
the method of Parisi-Sourlas. Since the effective action is second- 
order, the method of Ostrogradsky was used to quantize the system 
canonically. In the case of the harmonic oscillator, the relevant Green 
function can be computed exactly. In the general case, a perturbation 



O ■ expansion, involving time-dependent (memory) terms, can be defined. 
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In this paper we will show how to directly quantize open, dissipative 
systems that are subject to a random force. Previous efforts in this direc- 
tion were undertaken more than twenty years ago by Kostin and Yasue 
||. Kostin made use of a non-linear Schroedinger equation, which was later 
derived using Nelson's quantization |§ method by Yasue. Since then open 
systems have been quantized by "closing" the system first, which means as- 
suming a particular model of the environment and its interaction with the 
system {3|. Given the Hamiltonian for the system and environment, quanti- 
zation is straightforward. By integrating out the degrees of freedom of the 
environment, the quantum theory of the system is derived. 

The reason why a direct quantization of an open, dissipative system 
subject to a random force is not possible is, supposedly, the absence of a 
Langrangian or a Hamiltonian. However, as we will show in the next sec- 
tion, there is an effective Lagrangian appropriate to such a system. This 
Lagrangian can be derived using the method of Parisi and Sourlas 0. Start- 
ing with this Lagrangian, we quantize the system using the canonical path- 
integral framework. 

The equation of motion for an open, dissipative system subject to random 
force is 

d x dx dj^-J ^ (^i^J f 1 ) 

dt 2 dt dx 

where m, r\ and U(x) are, respectively, the mass, friction constant and po- 
tential of the particle. Because of the friction term, the LHS of equation ([!]) 
is not derivable from an action, i.e., there is no S a such that 

5S D d 2 x dx dU .g. 
5x dt 2 ^ dt dx 

For this reason, a direct quantization of the system using the canonical/path- 
integral method supposedly is not possible. 

But is there really no action that describes the system? Note that al- 
though there is no solution to equation the system is not described by 
(2) but by equation ([!]) plus the condition 

(f{t)f{t)) = i5(t-t). (3) 

Equation (H) implies that the distribution of / is gaussian, i.e., 

P[f]^eil dt ^. (4) 
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Following Parisi and Sourlas' arguments we find 

(x f (t)x f (t')) f = J (dx)x(t)x(t')e iJdt[ ^ ) +trln V^). (5) 

At the LHS, Xf{t) refers to the solution of equation (|l|) and the average 
over / makes use of the distribution given by equation (||). The RHS of 
equation (||) therefore is a classical two-point function of a particle in a 
random field. From this we read the effective action of an open, dissipative 
system subject to a random force as 

_ . . 1 r . ( d 2 x dx dU\ 2 . , ( 5 2 S \ 

So(x) = -Jdt^m— + r ] - + -j -itrln^—j. (6) 

This will be the starting point of the quantization process. 

Let us quantize a specific system. Take the case of the harmonic oscillator, 
with potential U(x) = \kx 2 . We can neglect the second term in equation 
(H) because ^Jf 1 is x independent. The resulting action is second-order and 
time-independent. The method of Ostrogradsky 0] is clearly applicable in 
this case. 

Define the coordinates and momenta 

<7i = x, (7) 
12 = x, (8) 

Pi = — ( ~tt~\ — T)[mx + rix + kx) — m[mg 2 +nx + kx], (9) 

ox dt \ox ) 

dL 

P2 = -7— = m \m'x + 772; + kx] . (10) 
ox 

The Hamiltonian is given by 

H = l^PiQi ~ L = -— p 2 q 2 P2Q1+P1Q2- (11) 

~[ 2 m z m m 

Note that the classical equations of motion in the Hamiltonian formalism 
given by 

dH , . 

Qi = 12 
dpi 

OH , , 

Pi = o > 13 

dqi 
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give the Euler-Lagrange equation, 







-±\ 


(-) 


dt 2 


v d'x J 


dt 


\dx 



, . i- a (14) 

Since one of the coordinates (q 2 ) is a velocity, it is a good idea to check the 
consistency of the Poisson brackets before proceeding with the quantization. 
The fundamental brackets are 

h,qj] = \Pi,Pj] = 0, i,i = 1,2, (15) 
[quPj] = Sij i,j = 1,2. (16) 

These brackets are consistent as shown by the fact that the equation of 
motion follows from them. Also the Jacobi identity 

[A, [B,C]] + [B,[C,A}} + [C, [A,B]} = 0, 

for any A, B and C function of qi and pj follows from the Poisson brackets. 

We can now quantize the system using the canonical/path- integral method. 
This results in the Hamiltonian path-integral given by 

{dqi ) {dq 2 ) (dpx ) (dp 2 ) e U dt [piqi +P2<?2 ~ H ^^ . 

The only minor complication in the derivation of the above expression is the 
operator-ordering of the p 2 q 2 term, which was resolved by symmetrization. 
Integrating out all terms except q\ = x yields 

PI= / \dx)e^ J ' ^(mx+vx+kxf _ ^ 



But this is just the path-integral given in equation (|5|), the only difference 
being the appearance of h in equation (jlTD- Thus, the stochastic dynamics 
of an open dissipative system is similar to its quantum dynamics. That this 
is true is supported by the work of Haba [0], who showed that quantum open 
systems are random classical dynamical systems. 

Since the effective action is quadratic, the two-point function is exact 
and given by the Green function of the operator m 2 -^ + (2mk — rf) + k 2 , 
which is 

I poo e %u(t-f) 

9{t ~ t,) = 2^ loo ^mV-^mfc-^ + P (18) 
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The contour evaluation of equation flHf) shows six distinct regimes: 





v 2 

m>k + -, 


(ii) 


m = k+ 2k> 


Hi) 


1 9 

— < m < — + k 
2k 2k 


(iv) 


rf , rf 
- 1 — - k < m < -!— 
2k 2k 


(v) 


m= 2k- k > 


(vi) 


^rf u 
m < k. 

2k 



For example, in the case of (Hi), the Green function is 



git - 



m (\ 



k 3 V4 



-a 2 /3 2 



a(t — t') + tan 1 



a 



where 



a 



P 



1 2mA; — rf 

2 + 4P ' 

1 2mA; — 7/ 2 

2 ' 



Next we consider potentials of the form 

1 



U(x) 



-kx 2 
2 



V{x) 



where V(x) is 0(x 3 ) or higher. In this case 

5 2 S d 2 d , d 2 V 

— j yj j j i 

5x 2 dt 2 dt dx 2 



(19) 

(20) 
(21) 

(22) 
(23) 



The second term of equation (|j) is time-dependent and represents memory 
of the particle's interaction with the environment as shown by the expansion 



trin 



s 2 s 

5x 2 



d 2 V 1 
dth(0) — (x(t))-- 



d 2 V d 2 V 
dtdt'h(t-t') — (x(t'))h(t'-t)- OI (x(t)) + . . . 



dx 2 



(24) 
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In equation fl2"4]), h{t — t') is the Green function of the operator m jm+V7fi, J t~k 
and is given by 

e iw(t-f) 



2tt 



The contour integration of equation 



-muj 2 + i^u; + k 
yields 



h{t - f) 



e m I 1 1 J 



r sm 



(4km— rj 2 )^ 



e m< Vsinh 



£ (4A;m - rj 2 ) 2 (t - £') 
£ (-4A;m + r7 2 )5 (i - f) 



(25) 

for 4A;m > ?? 2 , 

for4A;m = ?? 2 , 

ioi4km < rj 2 . 

(26) 

The propagators correspond to the underdamped, critically damped and 
overdamped cases. 

The derivation of the effective action shows that in the non-linear case, 
the system's memory of its interaction with the environment is given by 
explicit time-dependent terms. But in the case of the harmonic oscillator, 
the cumulative effect of the system's interaction with the environment is 
"hidden" in a propagator which is a Green function of a fourth-order operator 
777,2 SF ~^ (2mA; — rj 2 ) J^r + k 2 instead of the expected second-order operator 

To go around the problem of dealing with an explicitly time-dependent 
system, we use fermionic degrees of freedom to express the determinant in 
equation (El). The effective action we will deal with is 



S[x,<j>,<p] 



1 



d 2 x 



dx 



dV\ 



- / dt < m——r + r/— + kx + — — + m— — 



dt 2 



dt 

d 2 V" 
dx 2 



dx J 



dt dt 



(27) 



Again following the method of Ostrogradsky, the momenta are those given 
by equations fl9p,(|lOD and 



7I> 



7Tj 



dt 



ni- 



di 



(28) 
(29) 
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The Hamiltonian is 



H' = H+— 7I>7T^-— 7T^(7r^+77</)) + — (^+7/0)^ + — #71^ + [ k + ^ ) 0. 

m m m m \ ctar y 

(30) 

From this Hamiltonian, we construct the Hamiltonian path-integral, inte- 
grate the momenta and the coordinates 4> and <fi to arrive at the path-integral 
for x given by 

PI[x] = f (dx)e^ s ^f [x] (31) 



where 

n . , 1 f , / d 2 x dx , dV\ 2 . n ( d 2 d , d 2 F x 

s.//N = 2 y dt + ^ + fe + ^ J ^ + ^ + fc + 

(32) ' 

A perturbation expansion using the Green's functions g(t—t') and h(t—t') 
can be derived from equations ( PI) and fl32f) . 



In conclusion, let us put in context what was accomplished in this paper. 
To the extent that there exist classical, open, dissipative systems that are 
subject to a random force as given by equation ([I]), this paper can be used to 
directly quantize such systems. But equation ([!]) describes a "phenomeno- 
logical" system, which should be derivable from an underlying system and 
bath dynamics. 

The underlying dynamics is a matter of building models appropriate to 
particular systems. In general, different models will have different realizations 
of dissipative and random forces. For example, Caldeira and Leggett @] 
modelled the environment with a bath of different oscillators that have an 
unrealistic interaction with the particle (interaction does not fall-off with 
distance and not translation-invariant). In spite of these shortcomings the 
random and dissipative forces result from a quantum-statistical treatment of 
the system-bath dynamics when the oscillators have a particular frequency 
distribution. 

We now ask, should the "phenomenological" theory be quantized directly 
or should we start with the system + bath dynamics? Is the quantum theory 
of the phenomenological dynamics the same as that of the quantum theory 
of the system + bath dynamics then integrating out the bath degrees of 
freedom? In other words, are the two pathways (see figure) to the quantum 
theory of the system lead to the same result? In the case of the Caldeira- 
Leggett model, the two pathways give the same result as both yield the same 
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quantum ( X,yj) > (X w/ quantum friction 

and random force ) 



C A 



A 



classical ( X,y t ) 



->- 



( X w/ dissipation and 
random force ) 



Figure 1: Two pathways to quantum theory for an open, dissipative system 
subject to a random force (degree of freedom given by X). The bath degrees of 
freedom are specified by {yi). Starting from the underlying classical dynamics 
for (X, yi), the two pathways are given by (1) - steps (a) and (b); and (2) - 
steps (c) and (d) 

quantum Fokker-Planck equation. In the case of a slow degree of freedom 
interacting with a few, chaotic degrees of freedom the equivalence of the 
two pathways is not apparent. Berry and Robbins ||, using an expansion 
based on the small parameter e = ^ (ratio of the characteristic times of 
the slow and fast dynamics), showed that in the first-order there is classical 
dissipation but no quantum dissipation. This apparently violates the corre- 
spondence principle and possibly is an example of a system where the two 
pathways do not yield the same result. However, as we will show below, 
these observations are not correct. As for the apparent violation of the cor- 
respondence rule, a careful accounting of the two small parameters U and e 
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by Wilkinson |J showed that the classical dissipative force can be derived 
from the quantal Landau-Zener transitions. This means there is no violation 
of the correspondence principle. 

Similarly, it is not unreasonable to expect that the two pathways will 
yield the same quantum theory for the following reasons. Generally, step 
(a) involves a statistical treatment using the Liouville distribution function. 
Integrating out the bath degrees of freedom will give the reduced distribution 
function of X, which is used to compute the average of X, i.e., (X), and its 
dynamics. This should result in a classical dynamics which is dissipating 
and subject to a random force. Step (b), the direct quantization of an open, 
dissipative system, can be done either through the methods of Kostin, Yasue 
and this paper. 

Step (c), on the other hand, involves the density matrix, which satisfies 
the quantum version of the Liouville equation. Following this by step (d), 
which involves integrating out the bath degrees of freedom, results in the 
quantum theory of X. 

Since the classical limit of step (c) is step (a), then as long as we account 
properly the bath degrees of freedom in steps (b) and (d), the two pathways 
should yield the same result. The dominant configurations that give fluctua- 
tion and dissipation may not be the same for steps (b) and (d) as the Berry 
and Robbins example shows. In this example, classical dissipation is due to 
the transfer of energy from the slow degree of freedom to all the fast, chaotic 
degrees of freedom. In the quantum case, dissipation is a collective effect of 
avoided crossings described by the Landau-Zener theory. 

To substantiate the above conjecture, let us make use of the Berry and 
Robbins' example. Step (a) was carried out in the papers of Berry and Rob- 
bins and Jarzynski [[Kj]. Although the full Liouville distribution was not 
explicitly worked out (only the zeroth order distribution was given by Berry 
and Robbins in terms of the microcanonical ensemble while Jarzynski used an 
arbitrary function of the interaction hamiltonian subject to a consistency con- 
dition) the effective forces (Born-Oppenheimber, geometric magnetism and 
dissipative) acting on the slow degree of freedom were derived. Jarzynski's 
paper worked out the forces by deriving the Fokker-Planck equation satisfied 
by the reduced phase space distribution of the slow degrees of freedom. 

Steps (c) and (d), on the other hand, were worked out in detail by Berry 
and Robbins and by Wilkinson. Berry and Robbins' construction of the 
density matrix lends easy identification of the classical limit. Although their 
expression for the quantal (h finite) dissipative force in order e vanishes, 
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Auslaender and Fishman ITIfl showed that the classical value is derived if the 



limit ^ — > is taken first before t — > oo is taken. This step is equivalent to 
the reverse of step (b). 

In essence, the cited works show that steps (c), (d) and reverse of (b) 
is the same as step (a). It would be more desirable to show that following 
step (a) by step (b) gives the same quantum theory for the slow degree of 
freedom as steps (c) and (d). Since the effective dynamics of the slow degree 
of freedom does not involve a random force, the best way to carry out step (b) 
is by quantizing the classical Fokker-Planck dynamics derived by Jarzynski. 
Full equivalence of steps (c) and (d) to steps (a) and (b) is established if 
the quantum theory is the same as the Wigner distribution derived from the 
reduced density matrix. 

Finally, as a corollary to the above result, it must follow that if two baths 
characterized by the sets of degrees of freedom {yi} and {zi} yield exactly 
the same classical dynamics with dissipation and random force, then the 
quantum theory of the two system + bath dynamics followed by integrating 
out their respective bath degrees of freedom will also be the same. 
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